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KJQ' Abstract. The purpose of this paper is to further investigate the solution 

r^ ' space of self-similar spherically symmetric perfect-fluid models and gain deeper 

^^ ^ understanding of the physical aspects of these solutions. We achieve this by 

^^ , combining the state space description of the homothetic approach with the use 

of the physically interesting quantities arising in the comoving approach. We 
focus on three types of models. First, we consider models that are natural 
inhomogeneous generalizations of the Friedmann Universe; such models are 
asymptotically Friedmann in their past and evolve fluctuations in the energy 
density at later times. Second, we consider so-called quasi-static models. This 
class includes models that undergo self-similar gravitational collapse and is 
I- _ ' important for studying the formation of naked singularities. If naked singularities 

s«.i^ ' do form, they have profound implications for the predictability of general relativity 

Cn ' as a theory. Third, we consider a new class of asymptotically Minkowski self- 

Ca_^ ^ similar spacetimes, emphasizing that some of them are associated with the 

On ■ self-similar solutions associated with the critical behaviour observed in recent 

On ' gravitational collapse calculations. 
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1. Introduction 



Spherically symmetric self-similar (SSS) perfect-fluid solutions to Einstein's equations 
C^ ■ have attracted considerable attention over the last two decades. Although a nearly 

complete classification of such solutions has now been achieved, many of their features 
remain obscure. The aim of this paper is to provide a better understanding of the 
solution space and to unravel some of the physical properties of the solutions. 

There are a number of preferred geometric structures exhibited by the models, 
and these have led to several different approaches. Two of the most commonly 
used are the 'comoving' approach and the 'homothetic' approach. In the comoving 
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approach, pioneered by Cahill & Taub [Q, the coordinates are adapted to the fluid 
4-velocity vector, whereas in the homothetic approach, introduced by Bogoyavlensky 
and coworkers |0, p|, they are adapted to the homothetic vector. The main advantage 
of the homothetic approach as presented by GoUath et ah (P, ||; GNUl, GNU2) is 
that the state space is compact and regular. However, physicafly interesting quantities 
are more straightforward to obtain in the comoving approach, as shown by Carr & 
Coley (P, 0, ||; CC). It is therefore clear that the two approaches are complementary, 
and consequently both will be used in this paper. 

The outline of the paper is as follows. In section g we briefly summarize and 
compare the approaches of CC and GNU; the discussion here is primarily qualitative. 



with most of the mathematics being relegated to Appendix A . One new result 
presented in this section is the complete state space for the SSS models (see figure ^. 
In section ^ we focus on solutions of particular physical interest: the asymptotically 
Friedmann models; the so-called asymptotically quasi-static models, of which the 
naked-singularity solutions studied by Ori & Piran (g]; OP) are examples; and a 
class of 'asymptotically Minkowski' solutions. We also mention the 'critical' solution 
which has been discovered in gravitational collapse calculations, although we discuss 
this in more detail elsewhere [|lO[. A particularly interesting feature of SSS solutions is 
the way in which their characteristics change at certain critical equations of state and 
we list these equation-of-state dependent features in section ^. Section ^ contains some 
concluding remarks and indicates possible future research directions. The governing 



equations of both approaches are given in Appendix A , together with expressions for 
some physically interesting quantities. 

2. Geometric structures 

In this section, we summarize the approaches of CC and GNU and discuss their 
relationship. We will consider an energy-momentum tensor of perfect-fluid form 

Tf"" = liu'^u" + pigf"" - m'^m''), (1) 

where /x is the energy density, p is the pressure, u^ is the 4-velocity of the fluid and we 
choose units for which c = 1. The only barotropic equation of state compatible with 
the similarity ansatz is one of the form p = afi, where a is a constant [|l|. Causality 
requires — 1 < a < 1. We will conflne ourselves to positive pressure and the parameter 
interval < a < 1. 

2.1. The comoving approach 

In the spherically symmetric situation one can introduce a time coordinate t such 
that surfaces of constant t are orthogonal to fluid flow lines and comoving coordinates 
(r, 9, (f>) which are constant along each flow line. The metric can be written in the 
form 

ds^ ^ e^" dt^ - e^^ dr^ - R^ dn^ , (2) 

dn^ = d0^ + sin^ 9 dc/)"^ , 

where i/, A and R are functions of r and t. The equations have a first integral, m{r, t), 
which can be interpreted as the mass within comoving radius r at time t. Unless 
a = 0, this first integral decreases with increasing t because of the work done by the 
pressure. Spherically symmetric self-similar solutions can be put into a form in which 
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all dimensionless quantities, such as i/, X, S = R/r, /ii^ and M = m/R, are functions 
only of the dimensionless self-similar variable z = r/t M. Varying z for a given t 
specifies the spatial profile of various quantities. For a given value of r (i.e., for a 
given fluid element), it specifies their time evolution. 

Some of these quantities have a particularly straightforward physical interpreta- 
tion. Thus S is the scale factor, /i t^ specifies the density profile at time t and the mass 
function M is related to the divergence of the congruence of outgoing null geodesies. 
Values of z for which M = ^ correspond to a black hole or cosmological apparent 
horizon since this divergence is zero. Another important quantity is the function 

V{z)^e^-^z, (3) 

which represents the velocity of the spheres of constant z relative to the fluid. This 
should not be confused with the velocity of the fluid with respect to a Schwarzschild 
foliation, the 'radial 3- velocity', which we denote by Vr. Special significance is 
attached to values of z for which \V\ = ^/a and \V\ = 1. The first corresponds 
to a sonic point, the second to a black- hole event horizon or a cosmological particle 
horizon. 

It is convenient to introduce a dimensionless function x{z) defined by 

x{z) = {4Tr^r^)-"/^^+"\ (4) 

The conservation equations T^".^ — can then be integrated to give 

e-' = /3a;z2"/(i+«), e"^ = jx-'/"S^, (5) 

where /? and 7 are integration constants. The remaining field equations reduce to a 



set of ordinary differential equations in x and S (given explicitly in Appendix A.l). 
Thus, the governing equations in the comoving approach consist of an effectively 2- 
dimensional system of non-autonomous ordinary differential equations. However, it 
is also useful to regard the solutions as trajectories in the 3-dimensional {x, S, S) 
space (where a dot denotes zd/dz). Each trajectory is parametrized by the similarity 
variable z, so there is a 2-parameter family of spherically symmetric self-similar 
solutions for a given equation-of-state parameter a. Note that z may be either positive 
or negative but it always has the same sign as V. 

There is considerable arbitrariness in the selection of the axes that specify the 
solution space: although CC take them to be S, S and x, Foghzzo & Henriksen 



(|12|; FH) take them to be /j,, V and a variable linearly related to S. In either case, 
the functions have an obvious physical significance. The comoving approach thus 
has considerable intuitive appeal in that it affords immediate physical insights. For 
example, one can write down the metric explicitly and see immediately where the 
density goes infinite and the singularities occur. 

One disadvantage of the comoving approach is that the similarity variable z can 
be badly behaved. For example, sometimes finite values of z correspond to zero or 
infinite physical distances from the origin and sometimes zero values of z correspond 
to non-zero distances. Also many solutions span both negative and positive values 
of z, which means that z must jump form -|-cxd to —00. FH avoid this problem by 
using the coordinate ^ = 1/z but they then have a discontinuity at z = 0. Another 
disadvantage is that the comoving approach, at least with the variables used by CC 
and FH, makes rigorous proofs of general features of SSS solutions difficult. 

In {x, S, S) space the sonic condition \V\ = ^Ja specifies a 2-dimensional surface. 
Where a curve intersects this surface, the equations do not uniquely determine x (i.e.. 
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Figure 1. The different regions of the sonic line Q in terms of a V{z) diagram. 
Note that there exist two sonic lines with V = i\/a. The curves represent the 
subsonic parts (| V| < \AJ) of some solutions with a regular centre (to be discussed 
in more detail below). In addition, the static solution has been included. The 
thick curves correspond to the flat Friedmann solution and the static solution. 
Note that the orbits in one quadrant constitute the time-reverses of the orbits 
in the other quadrant. Thus, for example, the expanding Friedmann solution 
appears in the upper-right quadrant, while the contracting Friedmann solution is 
depicted in the lower-left quadrant. 



the pressure gradient), so there can be a number of different solutions passing through 
the same point. However, only integral curves which pass through a line Q on the 
sonic surface, the sonic line, are 'regular' in the sense that they can be extended 
beyond there. Thus the sonic surface 'filters out' the small subset of solutions that are 
physical. The behaviour near the sonic line can be treated as an eigenvalue problem 
P, 0|. It can be shown that the equations permit just two values of x at each point 
of Q and there will then be two corresponding values of V [|l^. If the values of V 
are complex, corresponding to a focal point, then the solution will still be unphysical. 
If they are real, at least one of the values of V must be positive. If both values of 
V are positive, corresponding to a nodal point, the smaller value is associated with 
a 1-parameter family of solutions, while the larger one is associated with an isolated 
solution. The eigenvector direction associated with the 1-parameter set is referred to 
as dominant or primary and the other direction as secondary. If one of the values of V 
is negative, corresponding to a saddle point, both values are associated with isolated 
solutions. 

On each side of a sonic point, i may have either of the two values. If one chooses 
different values for i, there will be a discontinuity in the density, pressure and velocity 
gradients. If one chooses the same value, there may still be a discontinuity in the higher 
derivatives of x. Only the isolated solution and a single member of the 1-parameter 
family of solutions are analytic (or at least C°°) and these are associated with the 
two eigenvector directions. One can show that the part of Q for which there is a 
1-parameter family of solutions corresponds to two ranges of values for z for each sign 
of z. One range (zi < \z\ < Z2) lies to the left of the Friedmann sonic point zp and 
includes the static sonic point z^; the other range (|z| > z^) includes the Friedmann 
sonic point. These features are illustrated in figure n^. When a = i, it happens that 
Z2 = zs and 2:3 = zp- 

CC have classified the p — a^ SSS solutions in terms of their asymptotic 
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behaviour. The key steps in their analysis are: (1) a complete analysis of the dust 
(a — 0) solutions, since this provides a qualitative understanding of some of the 
solutions with pressure in the supersonic {\V\ > ^/a) regime; (2) an elucidation of 
the link between the z > and z < solutions; (3) a proof that, at large and small 
values of | z | , all self-similar solutions must have an asymptotic form in which x and S 
have a power-law dependence on z; and (4) a demonstration that there are only three 
exact power-law solutions - the flat Friedmann model, a self-similar static model and 
a self-similar Kantowski-Sachs model. Prompted by results obtained by GNU, CC 
also find: (5) a family of power-law solutions which are 'asymptotically Minkowski' 
at large \z\; and (6) a family of solutions which are 'asymptotically Minkowski' at a 
finite value of z and which have a power-law dependence on In z. These new solutions 
only exist for a > 1/5; they were missed in CC's original analysis but they were then 
able to extend their work to include them. 

2.2. The homothetic approach 

When the homothetic Killing vector is non-null, the line element adapted to the 
homothetic vector field can be written in a diagonal form H 



ds -< 2T^^2^^2Trjj2^x)dT^ - dX^ ^ D2^{X)dn^ 



(6) 



where we have distinguished between when the homothetic vector field is spacelike 
{d/dX; upper line element) and when it is timelike {d/dT; lower line element). The 
unphysical spacetimes, ds^ and ds^ respectively, are hypersurface homogeneous [|I^ |T^ 
and the approach used by GNU exploits similarities with the equations governing 
spatially homogeneous models (see, e.g., jlq). 

By demanding that the energy density be non-negative, dominant quantities 
can be determined, and bounded dimensionless variables are obtained separately 
for the spatially and the timelike SSS regions. In each of the two cases, a suitable 
dimensionless independent variable is chosen, so that only one dependent variable has 
dimensions. Using these variables, Einstein's equations and the conservation equation 
T^'^-i, — result in a dynamical system in which the equation for the dimensionful 
variable is decoupled. The remaining equations form 4-dimensional autonomous 
dynamical systems, where the dynamical variables are related by a constraint (see 



Appendix A. 2). One of the variables is V, defined in the comoving context above. 
It is related to the tilt of the fluid flow with respect to the homothetic symmetry 
surfaces. 

Since GNU introduce four bounded variables, locally related by a constraint (for 



similar treatments, see |17| |18[), solutions are effectively treated as orbits in a compact 
3-dimensional state space. The fact that the state space can be compactified is of 
great advantage in that one can visualize all the solutions at a glance and this will 
be exploited throughout this paper. However, the physical interpretation may still 
be difficult. For example, parts of the boundary of the state space are not always 
themselves self-similar, as exemplified by the non-self-similar Kantowski-Sachs state 
space that appears as a boundary submanifold of the spatially SSS state space. 

The study of the self-similar spherically symmetric solutions now corresponds 
to studying the orbits in the state space of these dynamical systems. The future 
and past asymptotes of the orbits are important for understanding the corresponding 
models and are associated with equilibrium points in the state space. Indeed, there 
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Table 1. Interpretation of solutions asymptotic to the given equilibrium points. 



Label Interpretation 



C Solutions with a regular eentre or infinitely dispersed solutions. 

M, M 'Asymptotically Minkowski' solutions. 

K Non-isotropic singularity solutions. 

F Asymptotically Friedmann (isotropic singularity) solutions. 

T Exact static solution. 




Figure 2. The spatially SSS reduced state space for Q < i. For more details, 
see GNUl. 




Figure 3. The spatially SSS reduced state space for a > -p. Note that the 
equilibrium points ±M now have entered the spatially SSS state space. 
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Figure 4. The timelike SSS reduced state space for a < -r. The dashed curves 
indicate the boundaries of the sonic surfaces, located at V = ±^/a. The buUeted 
line in one of these surfaces is the sonic line Q. For more details, see GNU2. 



is a simple correspondence between these points and the asymptotic classification of 
CC, as discussed in sections 2^ and |[ Such points correspond to solutions with 
higher symmetry and always lie on the boundary of the state space. They often 
have a straightforward interpretation. The M points, for example, correspond to the 
Minkowski solution with a certain slicing. Here we will be interested in solutions 
asymptotic to the various equilibrium points, rather than the points themselves. 
Thus solutions asymptotic to the M points are 'asymptotically Minkowski' and those 
associated with F are asymptotic to the flat Friedmann solution. The interpretations 
of different asymptotes are summarized in table ll 

The labels used for the equilibrium points are sometimes accompanied by 
subscripts and/or superscripts referring to the values of some of the dynamical 
variables at that point. Note that models corresponding to an equilibrium point 
in a spatially homogeneous slicing, such as the flat Friedmann model, will correspond 
to orbits in state space in a slicing associated with the homothety. An orbit in state 
space describes a particular solution, each point on the orbit giving the solution on a 
particular homothetic slice (up to a scale factor, which is determined by the decoupled 
'scale' equation). 

Equilibrium points and flows along eigenvector directions for the spatially SSS 
reduced state space are depicted in figures || and g. There is a monotonic function for 
this system when \V~^\ > (GNUl). As there are no invariant submanifolds with 
V~^ = in the interior, all equilibrium points lie on the boundary of the state space. 
Equilibrium points and flows along eigenvector directions for the timelike SSS reduced 
state space are depicted in figures |4| and ^. As for the spatially SSS region, there is a 
monotonic function of the timelikc SSS state space when \V\ > unless \V\ — \fa H. 
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Figure 5. The timelike SSS reduced state space for a > -. 
equilibrium point M"*" now has left the timelike SSS state space. 



Note that the 



The sonic surfaces are located at \V\ — ^/a in the timeUke SSS regions, see 
figures H and 0. As discussed in section 2.1, solutions must cross the sonic surface on 
the sonic line Q in order to be regular. Such a sonic line is only present in one of 
the sonic surfaces. This, together with the fact that there is a monotonic function for 
\V\ > 0, implies that solutions entering the other half of the state space (the right half 
of figures H and g) arc unphysical (see GNU2). In addition, all orbits in the right half 
of figure ^(figure ^ pass through a sonic surface and end at the sink M"*" (M"*"), so 
that all these solutions have an irregular sonic point. It also turns out that solutions 
there have negative mass. 

In order to obtain a fully global picture of the solution space, the state spaces 
of the spatially and the tinielike SSS regions must be matched. This is done at 
equilibrium points for which \V\ = 1. The change of causality of the symmetry 
surfaces there corresponds to a Cauchy horizon in these variables. In figure ^ the 
spatially SSS state space has been matched with two timelike SSS state spaces along 
the lines of equilibrium points H. These lines of equilibrium points are artifacts of the 
homothetic approach and appear because the coordinates break down when \V\ = 1. 
Thus, the state space is 'pinched off' there, as seen in the figure. When matching, one 
must adjust the directions of the state space flows in the regions. For example, when 
following orbits from the timelike SSS region into the spatially SSS region via +II~ , the 
state space flow in the spatially SSS region will be opposite to that of flgures and 0. 
One of the matchings is done for V = +1. The sonic line of the corresponding timelike 
SSS state space will then be located in the sonic surface at y = -\-^/a. Consequently, 
due to the monotonic function discussed above, T^ < is an unphysical region in this 
state space. The other timelike SSS state space is attached where V ^ —1. In that 
state space, the sonic line will be located at F = — %/"; ^'^d V > is an unphysical 
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Figure 6. Matching the spatially SSS region with two timelike SSS regions. Note 
that the equilibrium points ±M^ are matched with M+ in each of the timelike 
SSS regions, as indicated by the dotted arrows. 



Table 2. Possible initial and final equilibrium points of orbits in the matched 
state space, together with the number of parameters necessary to describe the 
solutions near a certain equilibrium point. 

Sink/source Separatrix-generating One orbit 
Parameters 2 1 - 



a > 



±K"±, (M+) 



±F,CO 
±F, 0°, ±M 



T 
T 



region. 

The possible initial and final equilibrium points of orbits are the following (see 
table H): for a < | (see figures || and Q), the only equilibrium points that are sources 
or sinks are ±Kj_ in the spacelike SSS region and M+ in each timelike SSS region. Note 
that the points M"*" are located in the unphysical part of each timelike SSS region. 
In addition, -i-F and C" are separatrix-generating equilibrium points, i.e., points for 
which there is a 1-parameter set of orbits, spanning a surface in the interior of the 
state space. The static solution is represented both by the equilibrium point T and 
by an orbit through the interior of the timelike and spatially SSS state spaces. This 
orbit is the only interior orbit connected with T. When a > ^ the equilibrium points 
-l-M^ change stability from saddles and become sinks or sources. In addition, the 
separatrix-generating points ±M appear in the spatially SSS state space (see figure 
H). No new features appear in the left-hand parts of figures ^ and ^, i.e., the parts of 
the timelike SSS state space associated with physically interesting solutions. 



2.3. Comparison between the two approaches 

It should be stressed that some of the differences between the two approaches are 
not intrinsic but merely reflect the particular realizations used by CC and GNU. 



Self-similar spherically symmetric perfect-fluid models 10 

For example, one could in principle use compactified variables in the comoving 
approach and one could also express solutions in terms of the similarity variable in the 
homothetic approach. Indeed the variables used to represent the solutions in the two 
approaches must always be mathematically related in some way (except in spacetime 
regions where one or both sets of coordinates break down). However, since the 
coordinate descriptions of the spacetime are very different, physical comparison is non- 
trivial. Also, each approach suggests different 'natural' dependent and independent 
variables, so variables found in one approach can be difficult to express in the other. 

CC emphasize the behaviour of solutions in the asymptotic limit \z\ -^ oo, 
since the behaviour of solutions in this limit can only take one of a few simple 
forms (asymptotically Friedmann, asymptotically Kantowski-Sachs, and what they 
term asymptotically 'quasi-static'). However, the limit \z\ — > cxd does not usually 
play a crucial role in the GNU analysis but just corresponds to some region of 
their state space. In particular, there is no equilibrium point corresponding to CC's 
asymptotically quasi-static solutions since \z\ —> oo in this case corresponds to a 
surface within the state space. CC also discuss the behaviour of solutions in the limit 
z —^ and find that the solutions are either exactly static or asymptotically Friedmann 
at the origin. 

CC place considerable emphasis on the form of the velocity function V{z). This is 
useful if one wishes to identify event horizons, sonic points and physical singularities. 
(In the GNU analysis these conditions correspond respectively to the lines H, the 
surfaces |T^| = y/a and the equilibrium points K.) However, it must be emphasized that 
the V{z) representation alone yields an incomplete understanding of solutions since 
one is projecting them onto a particular 2-dimensional plane, so that many physically 
distinct solutions may be superposed. For example, only in the a = ^ (radiation) 
case does the region of the 3-dimensional solution space in which the mass is negative 
correspond to a well-defined region of V{z) space; otherwise it depends on the third 
axis. CC also emphasize the form of the scale factor S{z) and the mass function M{z). 
These quantities provide important physical insights into the solutions but one should 
bear in mind that they only represent different 2-dimensional projections of the full 
3-dimensional state space. 

In the homothetic approach, the variables are chosen for their mathematical 
properties. This makes it possible to obtain a compact state space, and thus a complete 
visual picture of the solution space. The homothetic formulation also provides rigorous 
proofs and enables one to use the insights gained from Bogoyavlensky's analysis 0. 
On the other hand, this approach is rather technical, and one can only indirectly gain 
physical insights. 

A complication in the diagonal homothetic approach is that the spacetime must 
be covered by several coordinate patches, since the coordinates break down when the 



homothetic Killing vector becomes null, |1^| = 1 (see section 2.2). This is physically 
equivalent to when the fluid 4-velocity becomes null. Since many SSS solutions cross 
\V\ = 1, one is forced to consider the relationship to other approaches if one wants 
to understand the global nature of each solution. Thus, in order to investigate when 
and exactly how orbits pass from one regime to the other (e.g., timelike to spacelike) 
in the homothetic approach, the dynamics are considered in the comoving approach. 
Going back and forth between the homothetic and comoving approaches leads to a 
complete understanding of the dynamics across the null surface and allows a check of 
the internal consistency of both approaches. This is another illustration of how the 
two approaches are complementary and how it is only by using the two approaches in 
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tandem that a full understanding is achieved. 

Variable and coordinate transformations can be used to relate the comoving 
and homothetic variables explicitly (see GNUl, GNU2 H, pi), although care must 
be taken in effecting these transformations, and it should be noted that there are 
some typographical errors in Appendix B of GNUl. An alternative and perhaps a 
more elegant procedure is to use expressions for dimensionless quantities, such as /ii^ 
and 2m/ R, in each approach. Together with the covariant expression for V, given by 
Cahill and Taub 0] (see equation A. IS), it is straightforward to use these to relate 
the variables of the two approaches. 

2.4- The Schwarzs child approach 

We should also mention the 'Schwarzschild' approach. This is useful when matching a 
self-similar interior region to an asymptotically flat (non-self-similar) exterior region. 
Schwarzschild coordinates are also useful if one wishes to solve the equations of motion 
for the null geodesies, as required in studying the global structure of the solutions. 
Consequently it was used by OP to study naked singularities and the cosmic censorship 
hypothesis in the context of SSS solutions. However, the disadvantage of Schwarzschild 
coordinates is that they involve non-physical singularities (see g], pp. 158-159). The 
transformations between comoving and Schwarzschild coordinates are given explicitly 
by OP. 

3. Solutions and their physical interpretation 

In this section, we discuss various types of solutions and their physical properties. The 
strategy will be to depict a given solution class in state space but - in order to extract 
physical features - to also plot the following physical quantities as functions of z: the 
scale factor S; the velocity function V\ the density profile ^t^; and the mass function 
2m /R (= 2M). Another quantity of physical interest is the asymptotic energy per 
unit mass E and this will be useful in parametrizing solutions. The expressions for 
these quantities are collected in Appendix A.3|. 



3.1. General solution structure 

Solutions can be classified by considering their global features, as well as their 
asymptotic properties. One possibility is to consider their differentiability. Another 
is to use the global causal properties of the homothetic Killing vector ry" and we now 
consider this. (1) There are solutions for which ry" is purely timelike; these develop 
a shock, since their orbits necessarily end at an irregular sonic point. (2) There are 
solutions for which 77" always is spacelike; these always have |y| > 1. (3) There are 
also solutions, such as the flat Friedmann solution, for which 77" changes causality once. 
(4) Solutions with two causality changes include, for example, the static solution and 
some asymptotically Friedmann solutions that recoUapse. (5) Solutions that undergo 
three changes of the causality include those that develop a naked singularity. For more 
examples, see table pi We will now discuss some specific examples in more detail. 

3.2. Asymptotically Friedmann solutions 

There are two distinct 1-parameter families of asymptotically Friedmann solutions as 
\z\ —> 00, one with z > and the other with z < 0. The parameter characterizing 
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Table 3. Classification of orbits. T and S stand for timelike and spacelike, 
respectively. 



Causality Example 



T Solutions with an irregular sonic point. 

S Recollapsing asymptotically Friedmann solutions, always having 

\v\>i. ^ 

Recollapsing asymptotically quasi-static solutions, always having 

TS Asymptotically Friedmann solutions that expand forever, including 

the flat Friedmann solution itself. 

Asymptotically quasi-static solutions that collapse to form a black 

hole. 
TST The static solution. 

STS Asymptotically Friedmann solutions that recollapse but have \V\ < 1 

for part of their evolution. 

Asymptotically quasi-static solutions that expand and recollapse, 

having \V\ < 1 for part of the evolution; these can also be of type 

STSTS. 
TSTS Asymptotically quasi-static solutions that collapse to form naked 

singularities. 



these solutions measures the underdensity or overdensity relative to the flat Friedmann 
solution and is also associated with the asymptotic energy E. In terms of the GNU 
state space, the asymptotically Friedmann solutions start at the F equilibrium points 
in the spatially SSS region (see figure |7|, where one of the two equivalent families of the 
asymptotically Friedmann solutions is depicted) . Physically interesting quantities are 
illustrated in figure |[ The z > solutions correspond to inhomogeneous models that 
start from an initial Big Bang singularity at z = cxd [t — 0) and then, as z decreases, 
either expand to infinity or recollapse. The z < solutions are just the time-reverse 
of these. Depending on the value of E, two qualitatively different types of solutions 
can be distinguished: 

1) For solutions that are sufficiently overdense with respect to the flat Friedmann 
solution (i.e., for E less than some critical negative value -Ecrit), \V\ reaches a minimum 
and then rises again to infinity as \z\ decreases, indicating the formation of a non- 
isotropic singularity at which 5 — > for finite a value of z. The solutions in this class 
are of type S or STS in the classification of section 3.1. Such solutions correspond 
to black holes growing at the same rate as the Universe [091 EG, O, pl|. They are 
represented by dotted curves in figures and ^. Providing the minimum of \V\ is 
below 1, there is a black- hole event horizon and a cosmological particle horizon where 
\V\ = 1. Otherwise the entire Universe is inside the black hole, although there is 
always an apparent horizon since CC show that the minimum of M is necessarily 
below 1/2. 

2) All solutions which are underdense or not sufficiently overdense {E > Edit) 
reach the sonic surface at \ V\ = y/a. The solutions in this class are of type TS in the 
classification of section 3T. Those which reach the sonic line Q, and in addition have 
zi < \z\ < Z2 or |z| > Z3 at the sonic point (see figure ||), may be attached to the origin 
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Figure 7. Asymptotically Friedmann solutions as orbits in state space. 
RecoUapsing solutions are represented by dotted curves. Also shown are solutions 
that reach the sonic surface (dashed curves), some of which can be continued 
to a dispersed state C. The continuations are exemplified by solutions with no 
oscillation in Vr (four curves) and one oscillation in Vn (one curve). The flat 
Friedmann solution is indicated by the thick dashed curve, 'x' marks an irregular 
sonic point. 
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Figure 8. Physical quantities for the asymptotically Friedmann solutions. The 
dash-dotted line in the V(z) diagram corresponds to the sonic surface. Other 
designations are given in the caption of figure M. The /jt^ diagram has been 
normalized so that the fiat Friedmann solution corresponds to \i^? = 1. The 
relation to the state space is as follows: the Friedmann point F corresponds to 
2 — » oo; the non-isotropic singularity (K) corresponds to 5 — > for finite values 
of z (dotted curves only); orbits reaching 2 = correspond to infinitely dispersed 
solutions (the point C; dashed curves only). 
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by subsonic (|V| < \/a) solutions. In figures andg, these asymptotically Friedmann 
solutions are represented by dashed curves. Subsonic continuations that have a regular 
centre have been included whenever possible. The solutions with a regular centre are 
also described by a single parameter and this is a measure of the density at the 
origin z — 0. These transonic solutions represent density fluctuations that grow at 
the same rate as the particle horizon [Ol . Numerical calculations indicate that these 
solutions can generally only be matched to the centre non-analytically. While there is 
a continuum of regular underdense solutions, regular overdense solutions only occur 
in narrow bands (with just one solution per band being analytic). The overdense 
solutions exhibit oscillations in the subsonic region, with the number of oscillations 
labeling the band. Solutions with larger number of oscillations form ever narrower 
bands within the one-oscillation band (in terms of the ranges of z at Q) . The existence 
of these subsonic bands was first pointed out by Bogoyavlensky and also studied 
by Ori & Piran [^, Carr & Yahil [Q and GNU2. The higher bands are all nearly 
static near the sonic point {z « zs), although they deviate from the static solution 
and approach the equilibrium point C as z ^ 0. 

3.3. Asymptotically quasi-static solutions 

As discussed by CC, there is exactly one self-similar static solution for each value of a, 
a 1-parameter family of solutions that are asymptotically static (in the sense that the 
radial 3- velocity Vr tends to zero as \z\ — > oo) and a 2-parameter family of solutions 
that are asymptotically 'quasi-static' (in the sense that S* — > but Vr is finite). One of 
these parameters can be taken to be the asymptotic energy E, while the other (denoted 
by D) can be related to the value of z at the Big Bang or Big Crunch singularity. It 
should be emphasized that the solutions in this class may only be close to the static 
solution for part of their evolution. For example, the class includes expanding models 
that finally recoUapse. This is clearly illustrated by some of the state-space orbits in 
figure ^ 

The key feature of these solutions is that they span both negative and positive 
values of z. Whereas the asymptotically Friedmann solutions are confined to z > or 
z < and are symmetric in z (changing the sign of z just reverses the time direction), 
the asymptotically quasi-static solutions are in general asymmetric in z and necessarily 
pass from z — — oo to z — +oo. This is because the Big Bang occurs at z = — 1/-D 
(corresponding to a non-zero value of t) in these solutions, so the limit \z\ — > oo 
has no particular physical significance. They can be interpreted as inhomogeneous 
cosmological models with an advanced Big Bang. Equivalently, for the time-reversed 
solutions, there is a Big Crunch singularity at z = -\-l/D. As illustrated in figures 
and on, there are two types of asymptotically quasi-static solutions: 

1) Expanding and recoUapsing solutions. Solutions with E less than some critical 
negative value Ecrit{D) expand from an initial singularity in the z < region and then 
recollapse to another singularity in the z > region. In terms of the classification 



of section |3J|, these are of type S, STS or STSTS. The value of z is -1/D at the 
initial singularity but depends on both E and D at the final one. These solutions are 
represented by dotted curves in figures ^ and |l^. As z decreases from —1/D, V rises 
from — oo, reaches a maximum below — -y/a but possibly above —1 and then tends to 
the quasi-static form as z ^ — oo. The solution then jumps to z = -|-oo and enters 
the z > regime. As z continues to decrease, V decreases to a minimum above ^/a 
and then tends to -f oo at the value of z corresponding to the recollapse singularity. 
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Figure 9. Asymptotically quasi-static solutions as orbits in state space. Dotted 
curves correspond to recoUapsing solutions, while dashed curves correspond to 
ever-collapsing solutions. In particular, the short-dashed curve corresponds to a 
solution having a naked singularity. Note that there also are solutions of the ever- 
expanding type. These go from the left-most K-point to the right-most C-point 
, in analogy with the ever-collapsing ones displayed here. The heavy full curve 
corresponds to the static solution. V and z are negative in the left half of the 
state space, and positive in the right half. 
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Figure 10. Physical quantities for asymptotically 
dash-dotted lines in the V{z) diagram correspond 
designations are given in the caption of figure bl The 
as follows: S ~> 0, fit^ —> oo, 2m/R — > oo correspond 
(K points); S ^ oo, V ^ 0, fit'^ finite, 2m/ R - 
dispersed solutions (C points; dashed curves only), 
the points K and when z jumps from — oo to oo, the 
a surface cutting the state space into two halves. 



quasi-static solutions. The 
to the sonic lines. Other 
relation to the state space is 
to non-isotropic singularities 
+ correspond to infinitely 
Note that |y| ^ oo both at 
latter case corresponding to 



Self-similar spherically symmetric perfect-fluid models 16 

If the niinimum of V is below 1, one necessarily has a black- hole event horizon and 
a cosmological particle horizon; this occurs if E exceeds some negative value E^:{D). 
The minimum of V will reach y/a when E reaches Ecrit{D) and this corresponds to 
the last recoUapsing solution. One also has the time-reverse of these solutions. 

The solutions with E > Ecrit{D) are of two types: 

2a) Ever-collapsing solutions. These start out from an infinitely dispersed state 
and describe the collapse of an inhomogeneous gas cloud to a non-isotropic singularity 



at z — -t-l/D (i.e., after t = 0). In terms of the classification of section 3.1, these are of 
type TS or TSTS. The dashed curves in figures g and |l^ correspond to these solutions. 
They start at the point C with V^ = at 2; = and then, as z decreases, reach a sonic 
point where V = —\fa. The solution is then attached to a supersonic asymptotically 
quasi-static solution at the sonic line. In this context, it should be noted that the 
introduction of the 'second' parameter D has relatively little effect on the form of the 
solutions in the subsonic regime. Indeed, one can show that all solutions apart from 
the exactly static solution must be asymptotic to the flat Friedmann solution at small 
\z\. In particular, the models can collapse from infinity (i.e., 5* — > 00 as z — > or 
t —f —00) only if E is positive or lies in discrete bands if negative. The supersonic 
solutions pass through a Cauchy horizon (where V — —1) before tending to the quasi- 
static form at z = —00 and jumping to z = +00. The continued evolution for z > 
resembles the evolution of case (1): as z further decreases, V first reaches a minimum 
and then diverges to infinity when it encounters the singularity at z = 1/D. The 
minimum will be below 1 if i? is less than some negative value E+{D), and in this 
case one necessarily has a naked singularity, as pointed out by OP. Particular examples 
of this are some of the general-relativistic Penston-Larson solutions (see figure 13 of 
OP). That D has different values for different solutions is clearly seen in the S(z) 



graph in figure 10 



2b) Ever-expanding solutions. These are the time reverse of the ever-collapsing 
solutions. They resemble case (1) in the z < regime but take a different form after 
they have passed into the z > regime. As z decreases from -l-oo, V (rather than 
reaching a minimum) decreases monotonically until it encounters a sonic point at 
V = y/a. If this sonic point is located in the physical ranges of the sonic line Q, the 
solution may be attached to the origin z = by a subsonic solution. The behaviour in 
the subsonic regime is equivalent to that discussed in case (2a). The ever-expanding 
solutions are not represented in figures |9| and |lO|. However, they can be obtained from 
the ever-collapsing ones (dashed curves) by letting z -^ — z and V -^ —Y. 

3.4- The asymptotically Minkowski solutions for a > -^ 

In this section we will discuss the 'asymptotically Minkowski' solutions, a class of self- 
similar solutions that were only discovered very recently and are of particular relevance 
to critical phenomena, see ^^ for further details. These were missed in CC's original 
asymptotic analysis, although they were able to extend their work after GNU had 
demonstrated the existence of such solutions. By combining the two approaches, we 
gain further physical insights into their significance. There are actually two families 
of asymptotically Minkowski solutions and both are only physical for a > i : 

A) The first family is described by two parameters and is associated with the 
equilibrium points M in state space. The solutions in this family have \V\ — *■ 1, 
S —^ 00 and // — » at some finite value z — z^:. In this limit, 

V/V = -{5a-l)/{l-a). (7) 
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Figure 11. Asymptotically Minkowski solutions of class A (dotted) and class 
B (dashed) as orbits in state space. Densely dotted and short-dashed curves 
correspond to regular solutions, while sparsely dotted and long-dashed curves 
correspond to singular solutions. 
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Figure 12. Physical quantities for asymptotically Minkowski solutions. The 
dash-dotted line in the V{z) diagram corresponds to the sonic line. Other 
designations are given in the caption of figure nil. The relation to the state 
space is as follows: S —> 0, \V\ —> oo, ^t^ — > oo, 2m/R —> oo corresponds to 
non-isotropic singularities (K points); S — > oo, V^ ^ 0, fit'^ finite, 2m/ R — > as 
2^0 correspond to infinitely dispersed solutions (C points); S — > oo, V finite, 
lit^ — > 0, 2m/ R — > with 2^0 correspond to the asymptotically Minkowski 
solutions (M, M points). 
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This family is illustrated by the dotted curves in figures O and O Although the 
value of z is finite at the M point, it should be noted out that the Schwarzschild radial 
distance {R = rS) is infinite. This is clearly seen in figure R2, where S goes to infinity 
for finite values of z. 

B) The second family is described by one parameter and is associated with the 
equilibrium points j-M in state space. The solutions in this family have |y| -^ K, > 1, 
S' — > cx) and /i — > as z ^ oo. The expression for T4 is 



a{a + 1) + y^a{a^ - a^ + 3a + 1) 

I/* - : • (») 

1 — a 

In figures |ll| and Q the dashed curves correspond to this family of solutions. 

Both families of solutions have 2m/ R — > as a solution approaches M or M. 
However, even though 2ml R —> as one approaches M, the mass m need not vanish. 
It is important to realize that these solutions are not asymptotically flat in the usual 
sense. Rather, they are perfect-fluid spacetimes for which the Minkowski geometry 
is obtained asymptotically along certain coordinate lines. This is reminiscent of the 
open Friedmann solutions, which asymptotically approaches the Milne model along 
certain timelines (see, e.g., ]l6[ 



As can be seen in figures nil and h3, both families contain solutions that can be 
connected either to the origin z = (via a sonic point) or to a non-isotropic singularity 
(for which S" — > at finite values zk). The former are of type TS, while the latter are 



of type S or STS, using the classification of section 3.1. For both families, the limit 
S ^> oo can be regarded as corresponding to an infinitely dispersed state, analogous to 
the late stage of an open Friedmann model, which is described by the Milne solution, 
see, e.g., Q. 

Note that there is a 1-parameter family of non-isotropic singularity solutions 
for each value of zk (i-e., one can take one of the two parameters that specify the 
solutions near K to be the value of the similarity variable at K). Most of these will 
be asymptotic to either a type A asymptotically Minkowski solution or a quasi-static 
solution. However, for sufficiently large values of |zk|, there will also be one type 
B asymptotically Minkowski solution and one asymptotically Friedmann solution, 
all connected with K. The limiting value of zk is the same in each case, reflecting 
the fact that the conditions at large values of \z\ have little influence on what 
happens near the singularity. This limiting value corresponds to the splitting of the 
asymptotically Friedmann and asymptotically M submanifolds: the orbits belonging 
to each separatrix surface can be described by one parameter. For some ranges of 
this parameter they are asymptotic to a K equilibrium point, corresponding to a non- 
isotropic singularity; for other ranges they have a sonic point. This is illustrated for the 
asymptotically Minkowski solutions in figure O, where dashed curves correspond to 
type B solutions. The analogous behaviour of the asymptotically Friedmann solutions 
can be seen in figure 0. 
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4. Examples of equation-of-state dependent features 

There are several features of the solutions that depend on the equation-of-state 
parameter a. In particular, the behaviour at the sonic line is greatly affected by the 
equation of state. Below, we list some of these equation-of-state dependent features. 



m a ~ 0.04: The boundary of the underdense subsonic band (see subsection 3.2) 
changes with increasing a from being the general-relativistic Penston-Larson 
solution to the degenerate node in the nodal region containing the flat Friedmann 
solution (i.e., z — z^). 

u a ~ 0.11: The upper boundary (as defined by OP) of the first overdense subsonic 
band changes with increasing a from being a secondary eigenvector direction in 
the nodal region containing the static solution to the degenerate node of that 
nodal region (i.e., z — 22)- 

• a = |: The equilibrium points M"*" leave each timelike SSS state space and enter 
the spatially SSS state space as the equihbrium points ±M. 

• a « 0.28: The critical solution, which corresponds to the lower boundary of the 
first overdense subsonic band, changes character with increasing a from a solution 
with a (generally irregular) second sonic point to a solution whose corresponding 
orbit ends at an M point in the spatially SSS state space. For the threshold value 
of a, the critical solution ends at an M point. 

• a = |: The degenerate nodes of the nodal regions of the sonic line coincide with, 
respectively, the flat Friedmann solution and the static solution (i.e., zp = Z3 and 
zs = 22)- 

The eigenvector direction at the sonic line of the flat Friedmann solution changes 
with increasing a from being a dominant direction to a secondary direction. For 
the static solution, the situation is reversed: the eigenvector direction going from 
being secondary to dominant. 

The boundary of the underdense subsonic band changes with increasing a from 
being the degenerate node in the nodal region containing the flat Friedmann 
solution to the flat Friedmann solution itself (i.e., z = zp). 

• a « 0.41: The eigenvector direction associated with the critical solution at the 
sonic line changes with increasing a from being an attractive eigenvector direction 
in the saddle region of the sonic line to a secondary eigenvector direction in the 
nodal region containing the static solution (i.e., zi < z < Z2). 

• a ~ 0.45: The overdense subsonic bands overlap, so that the band structure 
degenerates to only two bands: the underdense band and one overdense band. 

• a « 0.61: The zero in the radial 3-velocity Vu of the critical solution changes 
with increasing a from being subsonic to being supersonic (i.e., for this value of 
a, Vr = when \V\ — y/a). 

• a ~ 0.89: The eigenvector direction associated with the critical solution at the 
sonic line coincides with the degenerate node in the nodal region associated with 
the static solution (i.e., z = Z2), and changes with increasing a from being a 
secondary eigenvector direction to a dominant eigenvector direction. 
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5. Concluding remarks 

It is clear from the present work that we are close to a complete understanding 
of the self-similar spherically symmetric solutions. This paper has emphasized 
the advantages of a combined comoving and homothctic approach. This is par- 
ticularly clear in the investigations of the asymptotically quasi-static solutions and 
asymptotically Minkowski solutions. 

An important application is to study the critical solution of spherically symmetric 
perfect-fluid collapse. This has been done in a separate paper [Q, where it is shown 
that the critical solution is the unique self-similar solution which is analytic at the 
sonic point, has a regular centre, and contains one collapsing region surrounded by 
a dispersing exterior. As indicated in section ^ this solution is a member of the 
asymptotically quasi-static class for < a < 0.28, a member of the asymptotically 
Minkowski class of type B for a w 0.28, and a member of the asymptotically Minkowski 
class of type A for 0.28 < a < 1. 

Possible future investigations involve studying these models in more general 
contexts. One might ask if various features discussed here are structurally stable, 
i.e., how do other geometries and matter sources affect the behaviour? For example, 
one might further examine the role of the self-similar solutions in the complete space 
of spherically symmetric solutions. It also remains to consider the stiff fluid {a = 1) 
case, which is of relevance in the study of massless scalar fields. There are numerous 
bifurcations associated with the limit a — > 1 and this explains the numerical problems 
one encounters in this limit. However, we do not see any a priori analytic difficulties 
in studying these models using the methods presented in this paper. 
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Appendix A. Equations 

This appendix collects the basic equations in each approach, together with expressions 
for the physically interesting quantities that we have used. 

Appendix A.l. Com,oving equations 

The field equations (with G = c ~ 1) reduce to a set of ordinary differential equations 
in X and S: 

S + S+l-^^^--][S+il + a)S]^0, (A.l) 

\l + ai a X I 

^:^\ g4 , A^ ^2(l-a)/a^2(l-a)/(l+a) _ 

1 + aJ 13"^ S 

^^S^^(—^l\={l + a)x'^^-'^y'^, (A.2) 

X \ a I 
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M = 52^-(i+"'/" 



1 + (1 + «)| 



(A.3) 






(AA) 



where an overdot denotes z d/dz. At any point in the (x, S, S) space, for a fixed value 
of a, equations (A.3) and (A. 4) provide a constraint that gives the value of z; equation 
(A. 2) then gives the value of i unless \V\ = -Ja and equation ( |A.1| ) gives the value of 
S. Thus the equations generate a vector field (i, S', S") and this specifies an integral 
curve at each point of the 3-dimensional space, each curve representing one particular 

similarity solution. 

Where \V \ = \ /a, equations ( A.3 , A. 4) together with the expression for V{z) 
[see equation ( A 19| )] allow one to express S in terms of x and S, corresponding to a 
surface in (a;, S, S) space. Integral curves intersect \V\ = y/a in a physically reasonable 
manner only if 



f2aj'^ 






o4 , 2-^ (2-2a)/Q (2-2a)/(l+Q) 
^ +^25^ 



= (l + a)x(i-")/", (A.5) 

since otherwise the value of x and hence the pressure, density and velocity gradient 



diverge there. Since equation (A.5) corresponds to another 2-diniensional surface in 
{x, S, S) space, this will intersect the surface \V\ — %/oi on a line Q, the sonic line. 



Appendix A. 2. nomothetic equations 

Here, we briefly summarize the quantities and equations used in the diagonal 
homothetic approach. In the spatially self-similar case, the dependent variables 
(Qo: Q+j C*!, v) are defined in terms of the line element 



ds^ 



o2X 



Y2 



dr' 



(A.6) 



dx^ dn"^ ' 

together with Q+ = -D[/Di = Y'/Y + C{/Ci and v = 1/V, where ' = d/dr. Y is 
given by the decoupled equation 



Y' = -\q+ + Qo 



2Ql + 



1 



1 



:n. 



Y. 



The Friedmann equation defines fl„: 

The constraint takes the form 

G = (1 + a)wr2„ - 2 [1 + av^] Q+Ci 



0. 



(A.7) 

(A.8) 
(A.9) 
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Qo = - (1 - Oo') 



Q+'-C^^ + 



1 + Off 2 



i Lji 



Q+ — — <^oQ^ 



2(1 -Q+')- 



2n {1 + a)^„ 



1 + av^ 



1 {1 - a) + {3a + l)v' 



C{ =2Ci 



1 + av'^ 

Q+ + QoQ 



^ ^ni 



2 1 1 + a 
+ ^ 2 1 + a«2 



Qo^ri 



{(l + a)[2ago + (l + a)Q+]w 



(1 + a) [1 - ai;2] 

+ [a(3a + l)w2-(l-a)] Ci}. (A.IO) 

In the timeUke self-similar case, the dependent variables (S+, A, K^ F) are defined 
in terms of the line element 

3 



ds" = .^2T 






(A.ll) 



together with 1]+ = -D[/Di = 9' /O + A' / A, where ' = d/d^. 6* is given by the 
decoupled equation 



l + S+(l + E+)-i^-^^(i^f7. 



l + aV^ 



The Friedmann equation defines fit: 



Vlt^ 



{l-Y.\- A^ - K) 



The constraint takes the form 

G = (1 + a)V9.t - 2 [1 + ay2] Y.+A = 0. 
The reduced set of evolution equations is 



(A.12) 

(A.13) 
(A.14) 



S+' 



1 + aV'^ 



1 (3a+l) + (l-a)y2 



l + aV^ 



n. 



A' = 
K' ^2 

v 



1 + 2E+ + E+2 - A^ - -^^ TTTT^t 



l + aV^ 



A, 



tl-A^-^^^^^n, 

+ 1 + aV^ 



K, 



i-y2 



{l + a)[V^ -a 
+ [a(3a + l)-(l-a)y2] A}. 

Note that V was denoted by u in GNU2. 



{(1 + a) [2a+(l + a)E+]y 



(A.15) 



Self-similar spherically symmetric perfect-fluid models 
Appendix A. 3. Physical quantities 



23 



Here, we give the expressions for the different physically interesting quantities we have 
used. First we note that in the honiothetic approach, the similarity variable z ~ r/t 
is 



z ex S| 
ex SI 



;n(z;)e-^/'^-^"' 



(A.16) 

where equation (0) and the transformations between the comoving and homothetic 
approaches (see GNU1,2) have been used. 

The scale factor S = R/r is one of CC's dependent variables. In the homothetic 
approach, it is determined by 



S\Z) ex |;^|-(l-")/(l + ") 



IQ+I 
A 



(A.17) 



GNU use the velocity function (or its inverse) as one of their dependent variables. 
It can be covariantly expressed in terms of the fluid 4-velocity u" and the normal to 
the homothetic symmetry surfaces n" as follows m 



V 



^J{u°'naY - n^Ua 
In the comoving variables used here, V takes the form 

V = (/37)~'a;^'^"^/"^-'^^'-"^/^'+"^ 
The expressions for the density profile are 



(A.18) 
(A.19) 



fit^ 



J_2-22.-(l+")/" 

An 



S-'z- 



S-'z-' 



1 + av^ 1 — 



1 + aV^ K ■ 
The mass function is defined as M = m/ R where 



m{r,t) = ii? 



1 



-2i/ 



dR 
'dt 



-2A 



OR 
dr 



or equivalently 



m(r,t) ^Att iJLR^—-dr'. 
Jo or 



(A.20) 



(A.21) 



(A.22) 



In the comoving approach, it is given by either of equations (A. 3, A. 4), while in the 
homothetic approach it takes the forms 

1 + 2QoQ+ + Ql- Cl 



M 



2(1 -Qg) 
2K 



(A.23) 
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Finally, the asymptotic energy per unit mass E is the limit \z\ -^ oo of the 
dimensionless quantity £, whieh represents the total energy per unit mass of a 
comoving shell: 

£ 



12 



1 f [vjQo + Q+) + C, 

2\ K{l-V^) 



(A.24) 
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